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Abstract 



The hadronic vacuum polarization contribution to the muon (g — 2) value 
is calculated by considering a known dispersion integral which involves the 
R e + e -{s) ratio. The theoretical part stemming from the region below 1.8 
q^ \ GeV is the largest contribution in our approach, and is calculated by us- 

ing a contour integral involving the associated Adler function D(Q 2 ). In 
the resummations, we explicitly take into account the exactly known renor- 
malon singularity of the leading infrared renormalon in the usual and in the 
modified Borel transform of D(Q 2 ), and map further away from the origin 
the other renormalon singularities by employing judiciously chosen conformal 
CL|' transformations. The renormalon effect increases the predicted value of the 

O ' hadronic vacuum polarization contribution to the muon (g — 2), and therefore 

diminishes the difference between the recently measured and the SM/QCD- 
predicted value of (g — 2). It is also shown that the total QED correction to 
the hadronic vacuum polarization is very small, about 0.06 %. 
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I. INTRODUCTION 

The new precise measurement of the muon anomalous magnetic moment a^ = (g — 
2)/2 |l| allows for detailed testing of the standard model, and therefore for the possibility 
of looking into physics beyond the standard model as well. In fact, comparison of the 
experimental result with some theoretical calculations shows a 2.6a difference [JTJ . This has 
been suggested as the appearance of new physics.[] Since the advertised discrepancy comes 
from the calculation of Ref. || of the hadronic vacuum polarization contributions aX vp ) to 
Op,, we re-evaluate in the present paper the theoretical (pQCD+OPE) parts of this quantity. 
We use the same theoretical approach as in Ref. ||, which in turn is based on the approach 
of Ref. ||. However, in addition, we take into account the known renormalon structure of 
the Adler function. 

II. FORMALISM 



According to Ref. |10|| ) the hadronic vacuum polarization contribution to the muon 
anomalous magnetic moment a M = (g — 2)/2 appears in the following dispersion integral: 
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where K(s) is the QED kernel 

~2\ 



K(s) 




+ (l + xY 1 + 



x- 



ln(l+x) — x + 



x 



+ 



[1 + x] 



x lnx. 



x 



(2) 



Here, x = (1 -y M )/(l + y^) with y^ = (l-4ml/s) 1/2 . The largest part (about 92%) of a^ 
comes from the region with CMS energy y/s < ^/sq = 1.8 GeV. Following the approach of 
Ref. 0, applied in Ref. || to ay p '\ we rewrite the dispersion integral (JT|) (with s max = so) 
in the form 
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Here C\ is in principle an arbitrary constant, which however, according to the philosophy of 
Ref. H , may be chosen in such a way as to minimize the first ( "data" ) integral and maximize 
the second ("theory") integral. It is known that R e + e -(s) = 127rlmn(s + k:), with H(s) being 
the hadronic part of the (vector) photon vacuum polarization function which has no poles 



1 For an (incomplete) list of works investigating the possibility to explain and/or to find impli- 
cations of this difference within various frameworks of new physics, see: Refs. [Q (super symmetric 
models), || (non-supersymmetric extended Higgs sectors), § (composite models), [[| (models with 
extra dimensions), |§ (leptoquarks), and @ (other models). 



in the interval [0, 4m^); further, (1 — s/sq) also has no poles in that interval, in contrast to 
the function K(s). Therefore, the Cauchy theorem can be applied to the second ("theory") 
integral, with the path of Fig. [TJ. Carrying subsequently integration by parts, and using the 
identity D(Q 2 = — s) = — 12tt 2 s dll(s)/ds, leads to 
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where the associated vector Adler function D(Q 2 ) can be written in the following way:0 



D(Q 2 ) = N C ]T Q 
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Here, D can ,(Q 2 ) is the canonically normalized massless QCD part with dimension d = 0, 
whose power expansion in a (Q 2 ) = a ! f s (Q 2 )/'K is 
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with df = 1.6398 0, 4° } = 6.3710 p|, and df ] is estimated to be 4° } = 25 ±10 |3 
The renormalization scale in (0) is fi 2 = Q 2 , and the renormalization scheme is MS. The 



number of active quark flavors is nj = 3. The d = 2 contributions are [|l4j , |T5 
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where only the s quark contributes appreciably. The d = 4 contributions are those of the 
gluon condensate |b 
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and those proportional to m f [1 1 



2 This Adler function does not contain QED radiative corrections. The effects of them will be 
added later on. 
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The terms with dimension d > 6 do not contribute to the "theory" part in @ in the leading 
order renormalization group (RG) approximation. We note that the leading term in (^) is 
twice as large as that in || [their Eq. (9)]. We will see later that the terms ( |i"0"D give negligible 
contributions, but not the terms ©-©• In the quark condensate terms @ we can use the 
(approximate) leading PCAC relations m s (qq) ~ — f 2 {rn 2 K --rn 2 K ) and (m u +md)(qq) ~ — / 2 m 2 
(q = u,doTs; U = 0.0924 ±0.0003 GeV) 



«c) = E Q 2 f(m f ff) « -2.50 x 10- 4 GeV 4 

f=u,d,s 

(A<g)= J2 K/7) ^ -2.08 x 10- 3 GeV 4 . 
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Within the leading PCAC approach, the uncertainties in the numbers in fll~l)— (|12"D originate 
from the uncertainty of f 2 , and of the ratio e u = m u /m s = 0.029 ± 0.003 fli~6l|l5l . They are 
~1%, and affect insignificanly the results of the present paper. 

Insertion of the expressions ©-(PI]) m t° ^ ne con tour integral @ gives 
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Here, we used the complex momentum contour integrals 
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For RGE evolution of a MS (Q 2 ) we use the four-loop truncated perturbation series (TPS) 
[T7| of the MS beta function, with n/ = 3. In addition, for the RGE evolution of m s (Q 2 ) we 



use the MS four-loop TPS quark mass anomalous dimension ||18| . For example, the RGE 
evolution along the complex momentum contour gives 
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with 7i = 3.79167, 72 = 12.4202, and 73 = 44.263 [18 



III. EVALUATION 

We first use the input values as used in Ref. |§ 

(aGG) = (0.015 ± 0.020) GeV 4 , (18) 

afS(m 2 ) = 0.333 ± 0.017 (=> af^(M z 2 ) « 0.1201 ± 0.0020) , (19) 

m s (lGeV 2 ) = 0.20 ±0.07 GeV => m s (m 2 T ) « 0.151 ± 0.053 GeV . (20) 

Further, we use for D Cim .(Q 2 ) the NNLO TPS (i.e., with (ijj = 0) and with the renormal- 
ization scale /i 2 = Q 2 , i.e., the approach apparently used by [|J. The result for their input 
values, and for C\ = 0.007 GeV -2 , is then: 

10 11 x a^ p -\s < s ; th. part; NNLO L\ an .) = 4752 ± 108 , (21) 

in contrast to their value (4686.2 ± 113.2). The central values of these theory parts are thus 
higher by 1.4% than those given in ||. This percentage does not change at different values 
of the parameter C\, since the results are linearly proportional to C\. The uncertainty ±108 
in fl2T|) is obtained by adding in quadrature the uncertainty from a s (±61), from m s (±74), 
and from (aGG) (±49). The separate contributions to the central value 4752 in ([H]) are: 
4079 from the leading term; 749 from the (resummed) canonical part @ [ =^> (|14|)]; —87 
from the d = 2 strange mass term (|7|); 37 from the d = 4 gluon condensate term (§); —25 
from the d = 4 quark condensate terms (^); —1 from the d = 4 quark mass terms (0). 



We can, however, re-calculate the canonical part (14) by using methods which account 
for the renormalon structure of the Adler function. The latter structure has been shown 
fl3l , |T9"ll to have numerically significant effects in the hadronic r decay width ratio.0 Ref. (13 



uses ordinary Borel transforms, and Ref. [|I^] modified Borel transforms of -D can . In these two 
References, the exactly known leading infrared renormalon singularity of the Borel transform 
D can .(b), at the value of the Borel variable 6 = 2, has been fully taken into account as the 
pre-factor function of the form 1/(1 — b/2) 1+u and 1/(1 — b/2), respectively. This procedure 
allows us to describe the Borel transform more accurately in the most important region of 
the Borel integration, namely, the interval between the origin and the leading IR renormalon 
at b = 2. Furthermore, judiciously chosen conformal transformations b = b(w) have been 
employed there to map the singularities of other renormalons further away from the origin 
and thus to reduce their numerical significance. For details of the calculation procedures, 
we refer to these two References. 

If we re-calculate the caconical part (|T4|) by using the ordinary Borel transform method of 



Ref. |TS[], the result Q2"T| ) increases further to 4817 ± 116, as a consequence of the renormalon 
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In Refs. [13|,19 1, the massless QCD quantity r T was calculated (resummed). It is obtained from 
the hadronic r decay width ratio R T by subtracting from it the strangeness-changing contributions, 
factoring out the CKM-element |Vud| 2 an d the electroweak correction factors, and subtracting the 
quark mass (m u ^ ^ 0) contributions. The quantity r T is the same kind of contour integral in the 
complex momentum plane as the quantity ^4 ca n. of Eq. (|14|), but with the contour factor (1 + e 1 ^) 2 
replaced by (1 + e iy ) 3 (l — e iy ). Therefore, the analyses of Refs. |l3| , |l9| can be repeated, without 
any major changes, for the resummation of the quantity A can of Eq. ([lj). 



effects. In -D ca n. we took d 3 — 25, and for the renormalization scale (RScl) in D can we took 
the value which gives us the local insensitivity of the result with respect to RScl (/i 2 ~ 1.6Q 2 ), 
as argued in [O. 



If we apply to expression (|14D the method of resummation of Ref. |19| which employs 
modified Borel transforms, with the same input and with d 3 — 25 ± 10,0 we obtain a value 
very similar to the aforementioned one 



10 n x a™{s < s ; th. part; meth. [19]) = 4820 ± 120 . (22) 



The uncertainty ±120 in (|22| ) is obtained by adding in quadrature the uncertainty from 
a s (±74), from m s (±75), from (aGG) (±49), as well as from the resummation method 
uncertainty and the uncertainty of d 3 (±28). This calculation and the uncertainty estimates 
procedure are carried out in complete analogy with Ref. |19|| , to which we refer for details. 
The separate contributions to the central value 4820 in (|22|) are the same as those to (pi]), 
except that the contribution from the (resummed) canonical part ([]) is now 817 (before: 
749). 

In order to isolate the contribution of the renormalon structure included in the value 
(f2~2|), we should compare the latter central value with the one obtained by using for D can 
the N 3 LO TPS with d 3 = 25. The central value result in this case is 4772. The obtained 
renormalon structure effect is thus 48 units, or 1.0%. 

The input values (p^)-([20|), used by the authors of Ref. ||, and taken up until now in 
the present work, can be replaced by what we believe to be more updated values 



a™ s {m 2 T ) = 0.3254 ± 0.0124 (=> af S (M 2 ) « 0.1192 ± 0.0015) , (23) 

m s (m 2 T ) = 0.119 ±0.024 GeV . (24) 



The values (|23|) were obtained in |19| by a detailed analysis of the R T ratio, involving 
modified Borel transforms, and accounting for the renormalon structure of the associated 
Adler function via an explicit ansatz and with judiciously chosen conformal transformations. 



This result virtually agrees with the one obtained in the i? T -analysis of Ref. |13| where 
ordinary Borel transforms were used instead. The values (|23|) are shifted downwards and 
the uncertainties are reduced, in comparison to the values flH5|). The latter values are based 



largely on the ALEPH analysis of the R T decay |20| . The latter analysis did not account for 



the renormalon structure of the associated Adler function. The trend towards smaller values 
of a s and towards smaller uncertainties appears also in the analysis of the R T ratio of the 
authors of Ref. ||21|| , who accounted for the renormalon structure via a large-/?o resummation 



of the ordinary Borel transform and employed a resummation related to the effective charge 
(ECH) method - they obtained a s (m 2 ) = 0.330 ± 0.014. The values ([24]) for the strange 
quark mass, which are significantly lower than those in (|20|), were obtained in the recent 
analysis of Ref. Jl5[ . 



Further, ALEPH analysis |2(J of the r decays predicts the gluon condensate term to be 
consistent with zero 



4 We refer to [jDJ for a detailed discussion of this estimate of d 3 values. 



(aGG) = (0.001 ± 0.015) GeV 4 , (25) 



in contrast with the input flIB|). We will take, in addition to the (a s ,m s )-inputs 

either the input (|18 ) or (p5[) for the gluon condensate term. Small values of the gluon 



condensate close to the ALEPH values fl25"[) are also suggested in the formalism of Ref . [|22| . 



where the power-suppressed terms are obtained from the knowledge of the perturbation 
series of _D can (H) and of its infrared renormalon structure. 

Applying then again the resummation method of Ref. [|19|] to expression (|14]), we obtain 
the prediction 

10 11 x a^ p -\s < s ; th. part; meth. [19]) = 4823 ± 77 for (aGG) value Eq. flTJD (26) 

= 4789 ± 71 for (aGG) value Eq. (g5|) . (27) 



The uncertainties ±77, ±71 in (f26|)- (|27l) are obtained by adding in quadrature the uncer- 
tainty from a s (±49), from m s (±24), from (aGG) (±49, ±37), and from the resummation 
method uncertainty and d 3 -uncertainty (±25). The separate contributions to the central 
value 4823 in (|26|) are: 4079 from the leading term; 787 from the (resummed) canonical part 
(^); —54 from the d = 2 strange mass term (^); 37 [2 if using (p5|)1 from the d = 4 gluon 
condensate term (H); —25 from the d = A quark condensate terms @; —0.5 from the d = 4 
quark mass terms fllPl). 

We note that the uncertainties ±77, ±71 as obtained by us in (|26|)-(|27D are lower 
than the uncertainty ±113.2 for that quantity obtained by the authors of Ref. H (for 
C\ = 0.007 GeV -2 ). The main reason for this is that we used different, in our opinion 
more updated, values of a^ s (m 2 ) and m s (m^_) (|2"3"l)-(p4|) than the values (|T9l)-(|2"0"|) used by 
Ref. [§]. The uncertainties ±77, ±71 are about the same as the uncertainties obtained in 
Ref. for the "data" part, i.e., the first integral in (§, for C x = 0.000-0.007 GeV -2 (see 
Table. 1 of §). The authors of Ref. |] chose C x « 0.001-0.002 GeV -2 , i.e., their "theory" 
contribution was very small in comparison to their "data" contribution. The argument for 
the virtual exclusion of the theory from their considerations was that the uncertainties from 
their "theory" part were too high when C\ is appreciable. We believe that a different ap- 
proach, which emphasizes the "theory" part more than the "data" part, is legitimate as well, 
especially because the uncertainties from our analysis of the "theory" part are reasonable 
and comparable to the uncertainties of the "data" part even for large values of C\. The 
authors of Ref. || obtained the values of the "data" parts for C\ < 0.007 GeV -2 , by using 
the available e + e~ and r decay data. Therefore, we choose the largest C\ = 0.007 GeV -2 
listed in their Table 1. 

The "data" part given in their Table 1, for C\ = 0.007 GeV -2 , is 

10 11 x a^- p -\s < s ; data part) = 1622 ± 56 . (28) 

It corresponds to the first integral in Eq. (|j). 

Our evaluation of the theoretical part (|26])-([27l) of s < s contributions excludes any 
QED corrections to the hadronic vacuum polarization. The leading term of such QED 
corrections corresponds to the exchange of a virtual photon within the hadronic blob in 
Fig. H]. This diagram induces in R e + e - in the master formula (HD the real photon emissions 
and the virtual (plus soft photon emissions) corrections to the pure hadronic final states. 



Since the theory part in Eq. (f|) (2nd line) is dominated by the perturbative QCD (pQCD) 
Adler function, this QED correction can be easily implemented by replacing (1 + D CSuIi (Q 2 )) 
with the QED corrected one: 



in 
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(29) 



The QED corrections to the power-suppressed terms in (|5|) are negligibly small. Then, 
since the running of a ern (Q 2 ) along the integration contour can be safely ignored, the QED 
correction can be seen to shift the theory part by the factor: 
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Thus, the QED correction to the theory part of a 
It is small. 



(v.p) 



is approximately (6 x 10 4 ) x 4800 « 3. 



On the other hand, in the data part (28), the aforementioned class of QED contributions 
is, in principle, already included. This data part is taken from Ref. M, and corresponds to 
the first integral in Eq. (f|), for C\ — 0.007 GeV~ 2 , with R e + e -(s) there based largely on the 
experimental e + e~ and r-decay data. In this connection, we mention that there is some 
controversy on whether the specific photon decay channels via p meson [e + e~ —* (p) —>■ inc-/] 
are included in the data or should be added to them. The contribution of these channels to 
the data part (p8|), in the narrow width approximation, would be about (18 ± 4). However, 
the r-decay data already include this decay channel, while the situation with the e + e~ data 
in this respect is less clear. Since the authors of Ref. || included in their analysis both types 
of data, we presume that the additional possible contributions from the aforementioned 
photonic decay channels to the data part ( PB"D are significantly lower than (18 ± 4), and we 
will therefore ignore such contributions. 

We note that the inclusive, total QED correction to aL vp * s ver Y small when compared 
to the contribution from the real photon emissions. From the above QED correction to the 



theory part we can estimate the total QED correction to a[ t vp ' ) as a 
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x a e 
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4. This is much smaller than the contribution from the real photon 



7000 x (6 x io~ 4 ; 

emissions, which is according to [23[ about 90 units from the photonic decay channels of 
mesons, e.g. p, u, and 0. These two contributions may appear incompatible, but there is no 
real contradiction. It can be seen most easily in the narrow width approximation (n.w.a.) 
of the resonances. The n.w.a. formulas for the portions of R e + e - m (0) from the various 
decay channels of a resonance contain as factors the corresponding branching ratios, whose 
sum is always one. Thus, an increase in the photonic decay channel contributions must 
be compensated by the decrease in the pure hadronic channel contributions by exactly the 
same amount. The contributing QED correction factor in the Adler function in Eq. ( |29|) 
thus accounts for the QED corrections in the continuum only. Beyond the n.w.a. approxi- 
mation, however, the cancellation between the photonic channel contributions and the pure 
hadronic channel contributions cannot be exact, since the decay widths of the resonances 
are sizable, and so a small QED correction survives after the cancellation. This remaining 
QED correction, along with the QED corrections in the continuum, are those accounted for 
by the QED correction factor in the Adler function in Eq. (p9"l). 



We now include in our theory-part contributions (p6[)-(p7|) the aforementioned QED 
contributions (E9~|). 



10 11 x a^- p -\s < s ; th. part; meth. [19]) = 4826 ± 77 for (aGG) value Eq. flT3|) (31) 

= 4792 ± 71 for {aGG} value Eq. (H) . (32) 

This implies for the sum of the "theory" fl3~l~1)-(|32l) and the "data" ( |28| ) parts 

10 11 x a£- p -\s < s ) = 6448 ± 95 for (aGG) value Eq. (18) (33) 

= 6414 ± 90 for (aGG) value Eq. (|25[) . (34) 

The authors of Ref. |§ obtained for this quantity the value 6343 ±60, by almost excluding 
their "theory" part (complete exclusion of the "theory" part gave them 6350.5 ± 74). 

If we decreased parameter C\ from 0.007 GeV -2 toward zero, our central values in (|33])- 
(|3"4]) would go down toward the pure "data" value 6351, i.e., by ±97 and ±63 units, respec- 
tively. We will comment on this point in the next Section. 

If we add to the obtained quantities fl3"3"|)-(|3"4D the hadronic vacuum polarization parts 
from the region s > s as obtained in Ref. |8| ,f] we obtain 

10 11 x o {, v - p0 = 7029 ± 96 for (aGG) value Eq. (18) (35) 

= 6995 ± 91 for (aGG) value Eq. (0) . (36) 

The result obtained in the analysis of Ref. || is 

10 11 x a^- p - ] = 6924 ± 62 (D.H. [8]) . (37) 



While the bulk of the result (pFj ) of Ref. || was obtained by taking into account the data 
on e + e~ and the r decays, the bulk of the result fl3"5|)-(|3"6"D is obtained here by careful 
resummation of the contour integral of (^) where we account for the renormalon structure 
of the associated Adler function and for the d > terms. 



IV. COMPARISONS 

How do our results compare with recent results of others on oJ 4 v ' p '-'? 



5 They obtained 581 ± 15 for this contribution, using for the pQCD parts apparently the value 
( |lS| ) for a s . The authors of Ref. J23| obtained slightly higher values 584 ± 9, once we subtract 
from their values the (pQCD) contributions from (3 GeV 2 < s < sq) (so = 1.8 2 GeV 2 ). They used 
lower values for a s : o s (m 2 ) ~ 0.3088 ± 0.0245. This small discrepancy would apparently increase 
once we adjusted a s to the same value, say ([l9|), for the two approaches of Refs. |8],^|]). Then the 
prediction of Ref. @ would be about 592 ± 10. We decide to take the value of Ref. [§: 581 ± 15. 
We can reproduce their pQCD-parts with the simple TPS approach for R ee (s). When we use, 
instead, the approaches that account for the renormalon structure of R ee (s), and/or we use the 
value ( p3[ ) instead of ( |i"9| ) for a s (m^.), the values of these contributions change insignificantly. QED 
corrections are insignificant. 
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In Table | we show the values of a^- p ~> as predicted by others [p^ |2T|j23| ,|8[1 , along with our 
values.^] We see that our values are consistent with the values obtained by most of the other 
authors. We recall the fact that the bulk of our result is obtained by the pQCD+renormalon 
calculation, while the bulk of the results of the others was obtained by data integration. We 
observe that the consistency with the results of the others is even slightly stronger when we 
take the value of the gluon condensate of Eq. (pop, i.e., the small gluon condensate value 
obtained by the ALEPH analysis |20| and suggested also by the renormalon formalism of 
22fl . Several of the entries in the Table are based on inclusion of the r decay dataj] 

How do our results fl35|)-(|36|) compare with the experimental predictions for a„? This 



question remains somewhat unclear due to theoretical uncertainties of several higher or- 
der hadronic contributions, as argued by the authors of Ref. |23|]. The largest theoretical 
uncertainty is in the calculation of the hadronic light-by-light (1.1.) contributions. The 



chiral model (ch.m.) approaches would predict 10 
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x a 



(i.i) 



-86 ± 25 



the quark con- 



stituent model (q.c.m.) would predict +92 ± 20 [p3[] . The quark constituent model is valid 
for large values of the virtual photon momenta only, so we will give results for the chiral 
model unless otherwise stated. The QED corrections of the type of Fig. ^ have already 
been included in (^5])-(0). In addition, there are other QED radiative corrections where 
the photon propagator does not have both ends attached to the hadronic blob (the "rest") 
10 11 x a M (rad.corr.,rest) = —101 ± 6 [33]. These two radiative corrections (—86 ± 26) and 
(—101 ± 6) have to be added to the entries of Table | to obtain the total hadronic contri- 
bution a[ t hadr ^. In our case, we add them to the hadronic vacuum polarization contribution 



(we add the uncertainties in quadrature), and obtain 

10 11 x of adr ° = 6842 ± 100 for (aGG) value Eq. (18) 



6808 ± 95 for (aGG) value Eq. (g 



(38) 
(39) 



If we took, instead, the quark constituent model result for the light-by-light contributions, 
we would obtain 7014 ± 98 and 6980 ± 93, respectively. 



We can now add the results (p8|) (p9|) to the well known electroweak contribution ( |34,35 
and references therein) 



10 11 x fl/i (EW) = 116 584 858 ± 5 

and obtain the predictions 

10 11 x o M (predicted) = 116 591 700 ± 100 for (aGG) value Eq. (18) 
= 116 591 666 ± 95 for (aGG) value Eq. © 



(40) 



(41) 
(42) 



The actual experimental number |T| , when averaged with the older measurements |36| , is 



We did not include in the Table the results of some earlier analyses [28|. 



The entries of Table | are based on the use of the dispersion relation (|l]). The latter has been 



questioned in Ref. [29|, because there is no proof that the photon propagator is at least polynomially 
bounded. The inclusion of the r decay data for calculation of a^ ' has been questioned in Ref. |3(| , 
and the difficulties connected with this inclusion have been investigated in Ref. 
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10 n x ^(experiment, averaged) = 116 592 030 ± 152 . (43) 



The predictions (fiTD , (|42|) thus differ from the experimental result by 330 ± 182 (1.81a) 
and 364 ± 179 (2.03a), respectively, where a = 182, 179 is obtained by adding in quadrature 
the experimental (a 0X p. = 152) and the theoretical uncertainty (a t h. = 100, 95) .0 If we took, 
instead of the chiral model result, the quark constituent model result for the light-by-light 
contributions, the corresponding deviations would be 152±181 (0.84a) and 186±178 (1.04a), 
respectively. 

If we take only the newest experimental number [|]] 

10 11 x ^(experiment, new) = 116 592 020 ± 152 , (44) 

then the predictions (gl|)-(|42|) differ from it by 320 ± 182 (1.76a) and 354 ± 179 (1.98a), 
respectively. These deviations would be 0.78a and 0.99a, respectively, if we used the quark 
constituent model results for the light-by-light contributions. 



V. SUMMARY 

We obtain clear deviations of the theoretical results from the experimental ones when 
we use the chiral approach results for the light-by-light contributions. However, these 
deviations nonetheless are significantly smaller than the 430 ± 165 (2.6a) difference [|I[][] that 
has been suggested as the appearance of new physics. The deviations of our predictions 
from the new experimental values (|43|) , ([44]), range from 320 ± 182 (1.76a) to 364 ± 179 
(2.03a), depending mostly on the taken values of the gluon condensate. Our results are 
consistent with those of most of the other authors who used data integration for the bulk 
of their results. In contrast, the bulk of our results was obtained by pQCD+renormalon 
calculation. 

A major contribution to the aforementioned reduction of the deviations are our values 
for the hadronic vacuum polarization contributions (CLS1 and CLS2 in Table H), which are 
significantly higher than those of Davier and Hocker || (DH [|J in Table |). Our values 
originate from resummation of a major part of this contribution via a contour integral in 
the complex energy plane and accounting for the renormalon structure of the associated 
Adler function in the integrand. The accounting for the renormalon structure contributes 
about 50 units to 10 11 x al v ' p \ and this is larger than the expected future experimental 
uncertainties. 



8 The decrease of C\ to zero would contribute, in quadrature, additional 97 and 63 units to the 
uncertainties, respectively - as mentioned in the previous Section. In such a case, the predictions 
would deviate from the experimental result by 330 ± 206 and 364 ± 190, respectively. 

9 There, a ps 165 was obtained by adding in quadrature cr cxp . « 150 and a t h. ~ 67, where the 
value of cr t h. was taken from Ref . M . 
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-s = Q - plane 




FIG. 1. Application of the Cauchy theorem to the depicted integration contour in the Q 2 -plane leads 
to the contour integral in Eq. (0). 
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FIG. 2. QED correction to the hadronic blob. 



authors 



10 11 x a^ 



method 



ADH g§ 
BW H 



J [26] 



Nl 12 



N2 |3l 
TY1 (2 
TY2 jf 
DH | 
CLS1 
CLS2 



7011 ± 94 
7026 ± 160 
6974 ± 105 
7031 ± 77 
7011 ± 117 
7002 ± 65 
6982 ± 97 
6924 ± 62 
7029 ± 96 
6995 ± 91 



e + e~ + r data 

e + e~ data 

mostly e + e~ data 

e + e~ + r data 

e + e _ data 

e + e~ + r + spacel.F 7r (t) data 

e + e~ + spacel.F„-(t) data 

mostly e + e~ + r data 

mostly theory, and (|l8|) value 

mostly theory, and (p5[) value 



TABLE I. Comparison of predictions of the hadronic vacuum polarization contributions to a^ by 
various authors. 
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